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SUMMARY 

A surface integral is developed which yields solutions of the 
linearized partial differential equation for supersonic flow. These 
solutions satisfy boundary conditions arising in wing theory. Partic- 
ular applications of this general method are made, using acceleration 
potentials, to flat surfaces and to uniformly loaded lifting surfaces. 
Rectangular and trapezoidal plan forms are considered along with 
triangular forms adaptable to swept— forward and swept— back wings. 

The case of the triangular plan form in sideslip is also included. 
Emphasis is placed on the systematic application of the method to 
the lifting surfaces considered and on the possibility of further 
application. 


INTRODUCTION 

The study of supersonic lifting surfaces has been limited, 
until very recent times, to the work of relatively few theorists. 
The introduction of new propulsive systems has, however, increased 
the possible speed range of aircraft to such a degree that a great 
deal of emphasis has been placed on extending theoretical knowledge 
and techniques to cover controlled-flight velocities well beyond 
those in the subsonic regime. As a result of this added emphasis, 
within the past two or three years a number of important develop- 
ments and results have been reported which add materially to the 
main body of supersonic theory and included among these -publica- 
tions are several papers on wing theory for both the lifting and 
nonlifting cases. Most of these solutions, however, have been 
limited to the field of .conical flow, while a generalised theory 
which would fit any boundary condition likely to be met in problems 
of wing theory, has been lacking. 

In the present report advantage has been taken of the direct 
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analogy which exists between the linearized partial differential 
equation for supersonic flow in three dimensions and the two- 
dimensional wave equation of mathematical physics. As a result - 
of this correspondence, solutions which have been given for the 
wave equation are shown to be applicable to the type of boundary 
condition encountered in wing problems. The first section of the 
report is devoted to the development of the mathematical equation 
for the potential of the supersonic flow field. The application 
of this expression to a. number of examples in supersonic lifting- 
surface theory illustrates the usefulness of this method of solu- 
tion. In the first of these examples the loadings over the given 
plan forms are assumed to be' constant. The results obtained for 
such cases appear at first to be somewhat academic since undesirable 
twist and camber occur over portions of the resultant surfaces. 

From the constantly loaded surfaces, however, it is possible to 
develop surfaces having arbitrary load distributions . Imposing the 
condition that the final lifting surface shall be a flat plate leads 
to the solution of an integral equation in every case considered. 

The results obtained, for some of the plan forms considered, have 
been developed elsewhere but not with the unification of method 
attained here . New configurations are also included among the 
examples given. The methods shown are applicable to a large class 
of unsolved' problems of immediate interest. 
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LIST OF SYMBOIS 

(In order of their appearance) , 
velocity potential 
Cartesian coordinates 
local velocity of sound 

variable representing either the acceleration potential, 
the velocity potential, or any of the three perturba- 
tion velocity components 

free— stream Mach numbers 

transformed coordinates (See equation (3).) 
surface enclosing volume V 


V 


volume 
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) 


ni, ns, n 3 


P(X,Y,Z) 


r 

\ 


K 


direction cosines of normals to surface S 
direction cosines of conormal v to surface S 
surface on which boundary conditions are given 
point at which value of Q, i3 to he determined 
Mach forecone from point P(X,Y,Z) 
surface cutting axis of forecone T normally 
cylinder enclosing axis of forecone T 


u,v,w 


perturbation velocities in direction of X—, Y— , and Z-axes 
respectively 
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<Pu 

Po 

Pi 

Pu 

c o 

v 0 

pi 
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acceleration potential 

value of acceleration potential on upper side of . 
lifting surface 

value of acceleration potential on lower side of 
lifting surface 

density in the free stream 

static pressure on lower side of lifting surface, 

static pressure on upper side of lifting surface 

constant value of discontinuity in cp over uniformly 
loaded lifting surface 

free— stream velocity 

f —i T \ 

Mach angle of the free stream i p = sin- 

\ M 

span of wing 
chord of wing 

conical flow coordinate (See equations (27) and (30).) 
angle measured from X-axis 
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sn (u,k), j 
cn (u,k), \ 
d n (u,k) j 

F(u,k) 


K,K' 


E (u,k) 


E,E’ 


n (11,7) 



(3 tan 5 

angle of attack, radians 

pressure differential (p^— P u ) 

free— stream dynamic pressure (gP 0 V 0 2 ) 

load distribution 

lift of wing 

area of wing 

lift coefficient (L/qS 0 ) 
aspect ratio (b 2 /S 0 ) 

Jacobi’s elliptic functions of argument u and modulus k 

incomplete elliptic integral of first kind with argument 
u and modulus k 

complete e llipt ic integrals of first kind with modulus 
k and 1— k 2 , respectively 

incomplete elliptic integral of second kind with 
argument u and modulus k 

complete elliptic integrals Of second kind with modulus 
k and J 1— k 2 , respectively 

incomplete. .elliptic integral of third kind with argument 
u, parameter 7 , and modulus k 


Jacobi’s theta function 
i vertex angle of triangular wing 
A ■ angle of sideslip 

Hi, H 2 functions introduced in equations (105) and' ( 106 ) 


0 (u) 
A 
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THEORY 


L3.neariza.tion of Differential Equation for Compressible Flow 


The quasi— linear (i.e., linear in the derivates of highest order) 
differential equation for the velocity potential 0 in the case of 
compressible fluid flow in three dimensions, is expressible in the 
form 


/ 

$xx ^ 1 — 


®2c! 

a? 
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$y Z 
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$z 


— o 
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(1) 


where a represents the local velocity of sound in the medium and 
Cartesian coordinates are- used. Under the assumptions of small 
perturbation theory (references 1 and 2), this equation is modified 
so that it is linear in form and consequently more amenable to 
mathematical analysis. Denoting by the variable fi either the 
acceleration potential, the velocity potential, or any' of ‘the' three 
perturbation velocity components, the linearized ©repression for 
equation (l) is 


(l — M 2 ) Cixx + %y + ^zz " 0 


( 2 ) 


where M is the Mach number of the free stream and thus equal to 
the ratio of free— stream- velocity and the corresponding speed of 
sound. • , 


By means of the affine transformation . / 



i 

• Y = a/±( 1 — M 2 ) y j 


Z Ml - M 2 ) z 


( 3 ) 
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equation (2) can be put into standard forms. (The term affine refers 
to an arbitrary linear transformation of coordinates.) Thus, when 
M<1 the plus signs are chosen in the radicals of equation (3) and 
equation (2) becomes 

^XX + fi YY + fl ZZ = 0 00 


while for M >1 the minus signs are used and, as a consequence, 

%X ~ fi YY ” fi ZZ = ° 

For the case of subsonic flow (M<l) the linearized equation is 
thereby reduced to the well known Laplace equation in three- 
dimensions. Similarly, in supersonic flow (M >1) equation (2) is 
again reduced to classical type with the replacement of the space 
coordinate X by a time variable T to give the two-dimensional 
wave equation of mathematical physics. The linearization of the 
general differential equation for compressible fluid flow therefore 
makes available, in both subsonic and supersonic - studies, the results 
of the extensive work carried out in previous research on problems 
related to equations (4) and (5). 


Application of Green's Theorem to Linearized 
Compressible Flow Equation 

Methods of solution for partial differential equations of the 
type considered here may be classified inbo two principal categories: 
methods which express the solutions in terms of orthogonal functions 
and methods which are based on the use of Green's theorem. .Yol terra's 
solution, discussions of which may be found in references 3, 4, and 5, 
applies the latter approach to the two-dimensional wave equation and, 
as a consequence, his results may be adapted to the study of super- 
sonic flow and specific solutions of equation (5) . Let 


L(G) = Gjqj fiyy — 


The analytic form of Green's’ theorem for equation (5), relating 
a volume integral over the region Y to a surface integral over the 
surface. S enclosing V, may be written in the form 



OL(fi) — GL(o') ] dV = -/ f [cDn^ - n D n a l 


dS 
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where c, fl are any two functions- which , together with their first 
and second derivatives, are finite and single valued throughout the 
region considered, and 


ts n Bn 

D n n = n x ^ - n £ 


Bn Bn 
By " 113 Bz 


. where n x , n 2 , ns are direction cosines of normals to the surface S. 

The expression for D n Q is, of course, a directional derivative. 
The corresponding term appearing in Green's theorem for Laplace' 3 
differential equation (incompressible fluid flow) is precisely the 
directional derivative along the normal to the surface S. The • 
analogy "between the two terras prompts the introduction of the 
so-called conormal to S with direction cosines V x , v 2 , V 3 
defined as . 

v i = “ n .iy V 2 =t * 2> V 3 =r b 

The geometrical connection "between the normal and the conormal is 
indicated in figure 1; the angles between the lines and the Y— and 
Z-o.xos remain respectively equal, while the angles between the lines 
and the X-axis are supplementary. It follows, in particular, that 
if the surface S is the XY plane the two lines are coincident; 
if S is a cone with '•semi vertex, angle equal to 45° and axis parallel 
to the X-axis, then the conormal at any point lies along the surface S. 

It is now possible to write 


_D fi - v x — + V 2 — + v 3 s 2£L 

n Bx By Bz Bv 


( 6 ) 


and the surface— volume relation becomes 


r> n 


J 


I* j j, [ 01(a) - a i(a) ] dv = j I j 0 |fi. _ a 


Ba 

'SvJ 


dS 


(7) 


If fi and a are chosen so as to satisfy equation (5) 
throughout the region Y, then equation (7) reduces to the form 


n n 


r r ba P P Ba 

JJ a '8Z is= Jj 


(8) 
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The form of equation (8) la a direct analogue to results obtainable 
for functions satisfying Laplace's equation. (See, for instance, 
reference 6, page 46.) The use of the conormal produces this 
symbolic equivalence. 

Yolterra's Method for Threo-JXtinensional Wave Equation 

Consider now a surface r which, for the purposes of this report, 
may be thought of as being coincident with the XY plane and parallel 
to the air flow which is in the direction of the positive X-axi3. IVro 
such surfaces are represented by the darkened areas in figure 2. It 
is- desired to determine the value of at the point P(X,Y,Z) from 
a knowledge of the boundary conditions given on t. The solution to 
such a problem is immediately suggested by equation (8) since that 

equation requires only the knowledge of Q and along a surface 

enclosing a given volume, together with the knowledge of some 
particular solution c to the wave equation valid everywhere within 
the enclosed volume . Further, it is physically evident that in order 
that it contribute to the solution the enclosed volume should be 
within the forecone from the point P and also within the' envelopes 
of the aftercones from the foremost disturbances from P. Referring 
to figure 2(a), this would mean the volume bounded by the forecone T 
and the wedge X springing from the leading edge of r ; and in figure 
2(b), the volume . bounded by the forecone r and the aftercone X 
springing from the apex of the surface t, where in each case t 
represents the portion of the disturbance region affecting the value 
of Cl at P. Since for the boundary value problems involved the 
surface t remains in the XY plane, equation (8) must be applied to 
all three surfaces X, r and t. 

Since there is no way of determining fl and along, T the 

attempted solution must fail unless the particular solution a and 
its derivative with respect to the conormal vanish everywhere on F. 
But this is in fact the essential part of Volterra*3 method of 
solution. Thus the proper choice of a is 

. X-Xj. 

a = arc cosh — ■— ~ z 

J (Y-Yi) 2 + (Z— Zi) 2 

(This relation, incidentally, is the indefinite integral of the 
fundamental solution representing a 2 supersonic source in three 
dimensions [ (X—Xi ) 2 — (Y— Y x ) 2 — (Zr-Z x ) ] The value of a is 

equal to zero on the forecone p since the equation of this cone is 

(X-X 1 )2-(Y-Y 1 )2-(Z-Z 1 )2 = 0 
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and further, since the conormal is always directed along the forecone, 

^CT * 

— is the gradient of o along P and is also zero, 

dv 

Equation ( 8 ) now provides an equality for the distribution of 9 

and over X and r provided Cl and a satisfy eouation ( 5 ) 

dv ' ‘ 

throughout the enclosed volume mentioned. However, although a 
satisfies equation (5) everywhere in the enclosed volume opposite t 
from P— under the XY plane in figure 2 — along the line 
(Y— Yi) 2 + (Z— Z 1) 2 = 0 — above the XY' plane in. figure 2 — a is 
•infinite and does not satisfy equation (5) . If this line is excluded 
by means of a cylinder K of radius e, with axis lying along the 
line (Y— Y a ) 2 + (Z~Z\) 2 =0 , then equation (8) may be applied to the 
region outside If and yet within the space bounded by X, t, and r. 
In fact equation (8) can then be written 

nr 


i, U“ £ -a~)dS =0 ( 9 ) 

J j Tj+lf+X \ dv ov / 
where is the portion of t bounding the region of integration. 


IF P =v (Y— Yi) 2 -5- (Z-Zj.) 2 and cylindrical coordinates e, i|/: 
(X— X x ) are used., an element of area on the cylinder K is 
dS = — ed^d(X— Xj) , while 


and 


oo _ 80 _ 
bv ~ bTP. ' 


(X-Xi) 


so that 


e 7(X-Xi) 2 - 


lim 


lT f n bo „ ba \ 

. I Cl rr 9 ~ 'l' 

JJ xr \ uv ov y 


lim pp en(X-X 1 )dtd(X-X 1 ) 

e _5> 0 e Vli^T) 2 ”! 2 " 


lim 


1 


r, 


A X-Xi N 


i efl^arc cosh( ) a-^d.(X-Xi) 


- Cl tydXi -• P 


n *X 

2jt^ Cld, Y,Z)dt 


lim bCl f 


, III 

'j-JTZ g 0 c)v \ 


x 1 ~x'\ 


/ edtd(X— Xx) 


( 10 ) 


Appling this result to equation ( 9 ) gives 
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and. 



differentiating equation (11) with respect to X, yields 


n(x,Y,z) 


.-±LIT 

oX JJ t*j 


+ A 


5c 


a — ) dS 

ov 


( 11 ) 


( 12 ) 


Procedure for Lifting Surfaces and Symmetric Wings 

When the region considered is that hounded "by the surface t, 
r, and X 1 , the portion of X on the opposite side of t from the 
point P, then a is finite throughout the region and, as a direct 
consequence of equation (9) 


0 


i-4 ff — 

2jt oX Jw'r, + X»V ov ' Sv* 


dS 


(13) 


where ft’ is the value of the potential function on the side of t 
opposite P and v’ is in the opposite direction to v. 

Combining equations (12) and (13) gives 


fl(X,Y,Z) = LL ff 
2n <3X JJ Ti 



^£1 \ 
dv* J 


odS 


— — IT 

2rt 5XoUt-» 




5a 

5v 


dS 


ii rr ( a ^sV 

2 it oX JJ X V dv Sv / 


1 a rr / # 

2ir bx/A* V 




The integrations over r are now in a form which may he inter- 
preted directly in terms of known conditions over bodies with constant 

load or symmetrical section. The integration over X and X 1 can 
he handled by discussing the two cases shown in figure 2. If the 
leading edge is swept ahead of the Mach cone, as in figure 2(a), the 
value of Hi across the wedge springing from the leading edge is 
discontinuous, being zero immediately ahead of the wedge and a constant 
immediately behind it. It turns out, however, that in such cases the 
integration over X just cancels the integration over X’ provided 
only that the constants are equal and opposite above and below the 
wing. If the leading edge is swept behind the leading edge as in 
figure 2(b) the value of u, v and w are all zero on the cone 
(i.e,, continuous through the cone) and consequently, if Hi 


represents one of these perturbation velocities, HI, 2’, are 

all zero on X. Thus in either case there results the equation: 


(X,Y,Z) 


— — fP ( — + y ads - -i- — IT 

2 it Sx JJr x \ civ dV' 2 it 3 X'- / ~ /T i 




dS 

Sv 


(HO 


The counterpart of this equation for incompressible fluid flow is given 
on page oO of reference 6. 
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Under the particular conditions for which 


Sn = _ Ml 
5v 


over the surface t equation (ll) becomes 


n(z,Y,z) 



(15) 


( 16 ) 


The restrictions imposed in equation (ly) can he given physical 
significance after the functions Cl, Cl 1 and the surface t have heen 
given specific meanings. Consider first the case where t is a 
lifting surface. Obviously the normal induced velocity w is a 
continuous function across t. If Cl and fi* are velocity potentials 
associated with the lifting surface. 


w(I,Y,Z) -V(X,Y,Z) 

and equation (15) is satisfied. If q denotes acceleration potential 
or perturbation velocity u, it is necessary to show that on the 
lifting surface 


du _ du* 
5V SV 


This relation holds, however, for since w(X,Y,Z) = w’ (X, Y,Z) along 
r it follows that ' • 

3w _ idw* 

ax ~ ax* 


and from the condition of irrotationality it is possible to express 
the gradient of w in the X-direction as the gradient of u normal 
to the surface, that is, in the directions of V and v* . Equation 
( 16 ) is thus applicable directly to lifting surface theory in 
conjunction with either velocity or acceleration potentials. Applica- 
tion can also be made to the determination of pressure distribution 
over the surface of a symmetric airfoil at zero angle of attack. In 
this so-called nonlifting case the function Cl is set equal to the 
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induced velocity v, t is the plane of symmetry of the airfoil, 
and equation (l6) can be used to establish the boundary conditions 
provided equation (15) is satisfied. For this to be so bvfdv must 
equal -dw'/dv*. But conditions of symmetry give w(Z) = -*t' (— Z) 
from which the equality is seen to hold. 


Retransformation of Coordinates 


Since 


a 


= arc cosh 


X — Xj. 

V.(Y-Yi) 2 + (Z-Zi) 2 


direct substitution into equation (l6) gives 


Q (X, Y,Z) 


rr (fl-fl*) (a-X x ) (Z-Zx) dXydYi 

2«ZxJJt 1 [ (y_Y 1 ) 2 + (Z-Zi ) 2 ] V (X-X i ) 2 - (Y-Yi ) 2 - (Z-Z x ) 2 


This solution applies to equation (5) and, in order to relate 
problems to the linearized equation (2), it is necessary to use 
the transformation of equations (3). If the point' Xi, Yx, Zx 
transforms to the point x x , yx, Zx, it follows that 


l ^ TP (Q-fl*) (x-xx) (z-zx) dxi dyx 

^x Jj Ti [(y _ yi) 2 +(z _ Zl ) 2 ] VCx-xx)^^ (y-yxi^+Cz-zi)^] 


(IT) 


where 


f = M 2 - 1 
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APPLICATIONS 
General Remarks 

Applications in lifting surface theory may proceed along two 
possible lines depending upon the boundary conditions specified. In 
what is usually referred to as the direct problem the loading is 
given over the wing and the potential function of the flow field is 
calculated. From the potential function the shape of the aerodynamic 
surface supporting this load can be found relatively easily. The 
inverse problem concerns itself with the determination of the loading 
over a wing surface from a knowledge of the surface shape. In the 
following sections both of these problems will be considered. The 
direct problem will be discussed for various plan forms, the analysis 
proceeding directly from the expression for the potential function 
given in equation (17) . The detailed discussion of the direct 
problem is justified by its application to the inverse problem where 
the loading over flat plates with rectangular, trapezoidal, and 
triangular plan forms is determined. The mathematics of the inverse 
problem is less straight forward since the analysis involves the 
introduction of elemental lifting surfaces with constant loading 
and the solution of an integral equation, for each plan form. The 
cases considered are, however, of particular interest since the flat 
plate loading represents in thin airfoil theory the additional 
loading due to angle of attack which is associated with the 
particular plan form. 


Uniformly Loaded Lifting Surface 3 in Supersonic Flow 

Infinite span wing .— In order to determine the induced velocities 
on the surface of an infinite span, uniformly loaded, supersonic 
lifting surface by means of the methods derived in the preceding 
section it is convenient to set ft equal to the acceleration 
potential <p (reference 2) . The lifting surface is, in' this case, 
a surface of discontinuity for the function cp and corresponds to 
the surface r x in equation (17) . The discontinuity in the value 
of cp between the upper and lower surface is equal to 

(9 U -9,) - i (P, - P„) 

"0 


where 


density in the free stream 



Ik 


MCA TN No. 11H2 


3tatic pressure on lower surface 

P u static pressure on upper surface 

It follows that for the uniformly loaded wing in the plane 
z x =0 the discontinuity in the acceleration potential is a constant, 
say C Q . From equation (17) 


<pfoy,z) 


_1 ±rr 

2n c>xJJ 


C 0 (x-X!) z dX! dyx 

t (y-yi) 2+ z s W (3wci) 2 -p^[(y-yi)*+z s ] 


(18) 


A sketch of the airfoil plan form is given in figure 3 and two 
possible regions of integration are indicated. In all cases the 
integration with respect to y is performed between the limits at 
which the radical 

J (x-x x ) 2 — p 2 [ (y-y ! ) 2 + z 2 ] 


vanishes while the integration with respect to x depends upon the 
manner in which the forecone of the point P intersects the discon- 
tinuity surface. Denoting the chord length of the airfoil by c, 
the following relations are obtained: 

cp = 0 when x+0z< 0 

cp = ±iC Q when 0 <x+0z < c (19) 

cp = 0 when c < x+0z 

(When double signs are used, the upper sign refers always to the case 
where z >0 and the lower sign corresponds to z<0.) 

The value of the acceleration potential is thus seen to be zero 
at all points in space except for those points lying within the 
region between the wedges extending back from the leading and trailing 
edges of the airfoil. 

It is now possible to determine the induced velocities associated 
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with the acceleration- potential just obtained. Since, in linear 
perturbation theory (reference 2 ) 


Scp_„^u 9cp _ v dv .^P v dw 

5x " 0 ax' a y " 0 ax' al ' ° ax 


( 20 ) 


where u, v, w are respectively the x, y, z components of the 
perturbation velocities, it follows that 

. ‘ u = -i~ cp' 

V ° . ' ' 

y I zr <P ( x i>y* z ) dx i (21) 

ay J . v 0 

—03 


V 


a_ p* 




( X i>yxz) dx x 


The induced velocities for the infinite span airfoil result 
immediately from equations (19) and (21) , If the upper sign of a 
double sign . is again referred to the z >0 case, the results may be 
written in the form 



v = 0 


| 

\ for. 0 < x+’Pz, < c 
! 





( 22 ) 


Since the vertical induced velocities are constant, it follows that 
the supersonic airfoil of- infinite aspect ratio and uniform load 
distribution is a flat plate. The relations between this loading 
and angle of attack will be considered later. 

Lifting surface with rectangular plan form .— The complete 
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discussion of the supersonic lifting surface with uniform loading 
and rectangular plan form is lengthened considerably by the fact 
that in calculating the acceleration potential at the point P 
with coordinates x, y, z it is necessary to distinguish between 
several regions in space in which the point may be located. These 
regions arise from consideration of the manner in which the fore— . 
cone of. the point P cuts the surface of discontinuity. The value 
of cp can be found with approximately equal facility in each of 
these regions but, since this paper is concerned primarily' with 
effects on the surface of the airfoil, the solutions for pertinent 
regions only will be given here. 

Figure 1- shows the rectangular plan form LL’T r T together with 
the coordinate system to be used. The dimensions of the wing are 
chosen so that the Mach cones extending back from the leading edge 
will not intersect within the boundaries of the wing. This restric- 
tion, which is not necessary but merely simplifies the analysis, 
implies that if b is the span of the wing and c the chord length, 
then 


tan p = - - - <~ ( 23 ) 

Jw - 1 2c ' 

where p = arc sin i is the so-called Mach angle of the stream and 
equal to the semi vertex angles of the Mach cones. 

The loading over the rectangular plan form is to be uniform so 

the expression cp u -cp^ is set equal to C Q for — ib^yx^^b and 

O^Xx^c. The acceleration potential, expressed as a function of 
x, y, z, is thus obtainable from equation ( 17 ) and the limits of 
integration must be determined from the.position of P. From reasons 
of symmetry, only ,the portion of space from which y> 0 need be 
considered. Once the acceleration potential has been calculated, 
equations (21) may. be used to calculate induced velocities. The 
results of such calculations are given and the same convention for 
double signs is used. 

Region lx: Behind the leading— edge wedge, ahead of the trailing— 

edge wedge and bounded laterally by the inner sides -of .the two-tip 
Mach oones. The results in tMs region correspond to results 
obtained for the infinite span airfoil. 
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Thus 


<P = ±|C 0 
u = ±C 0 /2V 0 
v = 0 

V = - p 

2Yr> 


(24) 


Region I 2 : Within the Mach cone coning from the leading-edge 

tip hut outside the Mach cone from the trailing-edge tip. The 
integral of equation (17) must he broken into two parts. Denoting 
the integrand in the equation by the symbol I, 


1 r> T P^ 1 Ps^ p x 2 r> Y 2 1 

<P(X ’ 3 '- Z) “ - 5 W Q 4X t/ Tl “ yi t/r. tal X. Idyi J 


(25) 


where 


YX = 7 - | ^-Xi)2-p2 Z 2 X x = X - p </(y_ Jb)2 +Z 2 


Y 2 = y + | */ (x-x 1 )2-p2 z 2 x 2 = X^pz 


After integration of equation (25) and application of equations (21) 


q> 


= £° / + 1 


2 n 


arc tan 


X (y- £b) 


% V x 2 — P 2 [ (y— lt>) 2 +2« ! ^ 


u = 


a/.; 

2rtV 0 1 “ 2 


arc tan 


x (y- ib) 


z a/x 2 — p 2 [ (y— ib) 2 +z 2 ] 


l-i_ \o . _ o t ^ 


— z 


V = 


2«Vo (y- |b) 2 +z 2 


J x 2 — p2[ (y- ib) 2 +z 2 ] 


(26) 


arc xan 


P(y-i-b) 


w = (_ H + p 

2itV L 2 


+ t - - - — a/x?-P 2 [ (y- |b) £ +z 2 ] t 


V x 2 -P 2 [(y- ib) 2 +z 2 ] 


(y_ Ab)2 +Z 2 



18 


NACA TN Wo. 1412 


As a partial check on the expression for <p in equation (26) 
it can he seen that in the limit as z —> 0 the value of cp 
agrees with the result given in equation (24) on the wing while the 
value is zero off the wing. 

The values of vertical induced velocity in the plane z = 0 are 
of particular interest since from a knowledge of the distribution of 
w the surface shape and local angle of attack corresponding to the 
imposed load distribution can be determined. The expressions for w 
for uniform loading will be particularly useful later when the load 
distribution is modified in order to obtain airfoils with specified 
induced velocities. Introducing the notation 


(y- £b)g 

X 


(27) 


the following results are obtained for the area covered by the tip 
cone: 

(a) For x<c and — 1 <ti <+l 


w z=0 


£°L 

2nV 0 


it */ i— t] 2 

— + + arc sin p 

2 TI 


(28) 


For later reference it is desirable to have w z _q given in 
integral form. These expressions are 



Equations (28) and (29) indicate that the flow over the tip 
portion of the airfoil is of the type referred to as 'bonical flow." 
For this type of flow the values of induced downwash, aerodynamic 
loading, etc., are functions merely of the angle p. Busemann 
(reference 7), Stewart (reference 8) and Lagerstrom (reference 9) 
have developed analyses for certain plan forms which are postulated 
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on the existence of this type of solution and it appears that the 
theory can be shortened by such an approach. The chief value of the 
present method is most likely to be appreciated v, r hen no such assump- 
tions can be made. 

Tip of swept— forward lifting surface .- Consider the tip of a 
swept— forward supersonic lifting surface with uniform loading (fig. 5) 
the angle 5 Q between the leading edge and the X-axis and the angle 
63. between the trailing edge and the X-axis both being less than the 
free-stream Mach angle u. In carrying out an analysis it is neces- 
sary to distinguish between the type where the tip boundaries are 
behind the Mach cones and the type where the tip boundaries are 
ahead. The analyses of these two cases are of equivalent complexity, 
however, and can be handled with equal facility by the methods of 
this report. For all surfaces whose leading edges form an apex, 
only the case where the wing boundaries are behind the Mach cone will 
be considered. A Cartesian coordinate system is chosen as shown so 
that the origin lies at the vertex, the positive X-axis extending 
downstream, the Y-axis extending laterally, and the Z— axis being 
directed normal to the plane of the plan form and to the free— stream 
direction. The equations of the' sides of the lifting surface are 

y = 0 

y = - x tan S Q = - x 

and 

y = -(x-c) tan Ax = - (x-c) 

The calculation of cp(x,y, z) again must be divided into several 
cases depending upon the location of the point P(x,y,z) and again only 
solutions pertinent to the conditions on the lifting surface itself will 
be given. In the results listed below are included the general form of 
the integral for the acceleration potential cp(x,y,z) and the explicit 
expressions for ®(x,y, z) . The induced velocities, however, are given 
only in the plane z=0, as the integration to obtain~a general 
expression is difficult. The velocities in the z=0 plane, which are 
sufficient for the purpose of this investigation, can be obtained from 
a simpler integration since, for the integral involved, 

z ^ >0 yi(x,y,z) dx =y’l(x,y,o) dx 

This simplification was used in the analysis of most of the lifting 
surfaces investigated. As before, it is assumed that the discontinuity 
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in cp is equal to C 0 . Moreover, the expressions for v z _q are 
given in terms of the variables T and as where 

T 1 = and oi = 1 - - ( 30 ) 

In this manner the solution is shown to be conical in the region ahead 
of the trailing-tip Mach cone (fig. 5) • For points behind this Mach 
cone the flow is not conical but a function of both tj and oi. 


Region Ii: Inside the leading— tip Mach cone and ahead of the 

plane y + (0 o /p) x = 0 with y <0. From equation (17) 


<p(w) 


where 


/Xl 


nX g nO 

/ dXi/ Idyi +/ dx X/ 
L^o J Y X 


J *2 


pY 2 

Yi 




Yl 


^x 

p 


Y 2 = y - i V (x—Xi ) 2 — p 2 z 2 

p 


X, = 


(x+p0 o y)-P J (y+ ^ x) 2 +z 2 (l-0 o 2 ) 


1 - 0 O 2 


(3D 


x 2 = x - p Jy 2 +z s 


After integration 


2tt 


— arc tan 


ry 


z V x 2 — p 2 (y 2 +z 2 ) 
In this region, — 1< Tj <—6 q . 


+ arc tan 


xy+x 


2 


T 


- pz 0 


2„ -1 


z V x 2 — p 2 (y 2 +z 2 ) 


w 


z=0 


„ _ W> f n l 

li ii 2 (ii + e 0 ) 


2rtVr 


dTii 


(32) 


(33) 


C 0 P T 1— 1 ! 2 1 -1 / 1— 0^ 2 

arc cosh — + 

2rtV n Ti 0« T 0, 


-2- arc cosh 


-(i + e 0 n) 

(0 


+ 0o' l O ~ 

n +T l ) - 



I 


. NACA TN No. lhl2 


21 


Region I 2 : Inside the leading-tip Mach cone, outside the 

trailing— tip Mach cone and he ween the planes y = 0 and 
y + (0 O /|3) x = 0. Using the notation of equation (32) for 
etc., the acceleration potential is 

/ n Xi p o r X 2 p° p% 3 r> ■‘•Q 

l o a Vy, Iiy 4 1 tal 4 Ily 4 2 “4 I4yi . 


2 it dx 


X 1j Yi, 

,Y a 

“4 


where 


m 


Y 3 = y + j J (x-x 1 ) 2 -p 2 z 2 
X 3 = x+pz 


Integration yields 


x*0. 


<P = 


2n 


— arc tan 


*y 


xy 


— Pz 2 8 


+ arc tan 


z (y^+z 2 ) 


Zv />^-e 2 (:4+z 2 ) 


-] 


( 35 ) 


For the region considered — G c < T|<0 and in this case 



^ 1-3 1 2 


ni 2 (Tli + 0 o ) 




( 36 ) 


C 0 P 

2*V 0 


f ^1-T] 2 1 

l n A 


arc 


-i Vi-e 0 2 

cosh — t — 

1 K 


arc cosh 


(i+e 0 r)) 1 

(Q 0 +n) J 


Region I 3 : Inside the leading— tip Mach cone, outside the trailing— 

tip Mach cone and for values of y greater than zero. It follows that 


cp(x 


* r nXi p o pX 2 ro 

,y,s) = r~ / dx i/ Id Yi +/ dx i/ I d Xi 

2irdxWo d Yi '>X x j Y 2 j 


( 37 ) 
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The integration yields 


q> = 


Co 

2jt 


— arc tan 


xy- 


z Vx 2 - 3 2 (y 2 +z 2 ) 


+ arc tan 


y?e a 2 

xy + — p 3 Z 0 o 

Z fj x2_p2 (y 2 +z2 ) 


(38) 


In this region, 0 <rj <1 


v z =0 


-c^ 4 

2«V 0 J x ‘ t 1i 2(t 1i + 0 o ) 


_ CoP f 

2rtV r ' 


fv 1— T] 

L n 


— arc co3h — + 
0n "H 


/.^° g arc cosh iilMi 
e 0 ( 0 o +T i) 


(39) 


Region : Inside both tip Mach cones and ahead of the plane 

y - — ( 81 /P) (x— c) for y <0. The solution in this region is greatly 
simplified by use of the fact that for linear differential equations 
any algebraic sum of solutions will be another solution of the equa- 
tion. Since the differential equation for the acceleration potential 
is linear, this property can be applied to obtain a solution for this 
section of the airfoil by subtracting from the expressions given for 
region I 2 the expressions given for region I x except that in the 
latter region the variables are changed such that 5 X replaces & Q 
and x— c replaces x. From this operation it follows: 


cp = 

2n 


arc tan 


«xy 


+ arc tan 


0n 

xy + x 2 -30 o z 2 

+ arc tan E , , 

Z V X 2 -P 2 (y 2 +Z 2 ) Z - 7 x 2 - "pS (y2+ Z 2) 

(x-c) y (x-c)y+ ^(x-c) 2 -36;iZ 2 

~ arc tan — ■ ■ 5 — 

z x— c) 2 — 3 2 (y+z 2 ) zV (x-c) 2 -3 2 (y 2 +z 2 ) 


(40) 


w z=0 = 


Co3 


_ ^o 


2irV 0 


■— T : 2 1 . -1 a/ 1 - 0 O 2 (l+ 0 o Tl) 

arc cosh — + arc cosh 


So ^ e c 


(0 O +T l) 


— v ' /(ug ~ 7 1 2 + 1 


arc cosh 


f _ Co'S _ 1—0 1 ^ 


V 11 


arc cosh 


01 


V 0iP rfr l 


(41) 
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Lifting surface with trapezoidal -plan form .— The linear 
property of the differential equation may he used to advantage in 
determining the flow about a trapezoidal lifting surface with 
uniform lift distribution, since the boundary conditions within the 
plan form of the airfoil are obviously satisfied when the accelera- 
tion potential for a triangular tip i3 subtracted from the potential 
for the rectangular surface. 

Suppose (fig. 6) the angle of rake of the trapezoid is 
and that 5 0 is less than the Mach angle u . The acceleration 
potential will be identical, over the central portion of the surface, 
to that for the lifting surface of infinite aspect ratio. Over the 
parts of the surface which are blanketed by the tip Mach cones the 
flow will, however, be modified. Because of symmetry the determine— 
tion 01 this modification need only be carried out on one side of 
the figure. 

If the coordinate axes are chosen as shown in figure 6, the 
lateral boundary of the lifting surface is 

y = — x tan 8 0 - — ^25 

0 


It has been shown that both the rectangular plan form and the 
triangular plan form experience conical type flow over the region 
within the tip Mach cones. Thus, the variable rj defined in 
equation (30) may be used. 

Region lx : Inside the Mach cone originating at the leading- 

edge tip and ahead of the plane y = - |& x (i.e., on the surface 
of the wing), - 1< ’]<— 9 Q : 
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Region I 2 : Inside the Mach cone originating at the leading-edge 

tip, outside the Mach cone originating at the trailing— edge tip, and 


between the planes y = 0 and y = - 


- 6 n <r\< 0: 


C 0 (3 f x/l-Pi 2 cL'Oi 

z -° 2itV 0 \ J_ 1 Til 3i + 


\ 

T ll +0 n j 


c 0 P 


2rtY, 


+ arc tan 

. 2 


T! 1 

. — + — arc cosh l 

Vl— Tl 2 V T) / 


/i\ 7 i-e Q s 

)sh 


arc cosh 


( 43 ) 

1+0^ ‘ 


0 o + 1 > 


Swept-back lifting surface .— As another example of the way in 

which the linearity of the differential equation may be utilized 
to obtain further solutions, the induced vertical velocities for a 
swept— back wing will be determined for the case in which the leading 
and trailing edges lie behind their respective Mach cones (fig. 7 ). 
The boundaries of the plan form shall be given by the equations 




x 



y = 


o i 
(3 


(x-c) 


The flow will he conical ahead of the trailing— edge Mach cone 
where the induced velocities can he expressed in terms of the 
variable tj. Behind the trailing— edge Mach cone the flow will not he 

conical hut will he expressible in terms of the variables and 

co = 1 - 

Consider first the region of conical flow. In order to determine 
w z _o for a given value of p it is possible to consider separately 

the induced effects reduced by each half of the surface. But in 
this region ahead of the trailing— edge Mach cone, the induced 
velocities arising from one half of the surface are given by the 
formula for a similar region on the swept— forward surface. For 
reasons of symmetry the results for the entire swept— back lifting 
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surface need only be given for values of n within the limits 
- 1 <Tj <0. 


In the region where the flow is not conical the solution will 
he built up of a combination of solutions obtained from the regions 
of conical flair. . 


Eegion : Inside the leading— edge Mach cone, outside the 

6 

trailing— edge Mach cone and ahead of the plane y = — ~ x. From 
equation ( 39 ) the value of v^q resulting from the induced effects 


of the left-hand portion of the triangle can he determined and from 
equation ( 33 ) the induced effect produced hy the right-hand part of 
the triangle is given provided the sign of rj is changed. Combining 
these two results , it follows that 




z=0 : 




r 

| 


&tV r 


arc 


. -1 
cosh — 


vi-e 0 2 f 
+ 1 arc 

0~ V 


cosh 


1+M ^ 

+ arc cosh 

0 O +T 1 


1-0 o t] 


= c oP ( 
2*V 0 \ 


2 _ 

0n 


arc 


cosh — 




arc 


coah 


2— 0 q 2 -t1 2 '\ 

0Q 2- T l 2 ) 


m 


Eegion I 2 : Inside the leading— edge Mach cone, outside the 

n 

trailing— edge Mach cone, and between the planes y = — _£ x and y = 0-. 

3 

In the same manner as in the preceding case equations (39) and ( 36 ) 
can be combined to give this solution. Thus : 


w 


C 0 P f 


2=0 2itV 0 l 0 O 


-1 Jl-Q* 

nVi - 4 - _ _ __ u 

r l 0Q 


— - 1 — 0 qI) ( l+fl-T) 

arc cosh — + — [arc cosh + arc cosh— \ 0 


C Q P f-2 -1 , v 

= — — I — arc cosh — + 

2*v 0 ve 0 T! 0, 


k-°o 2 


arc cosh 


2-e 2 -t)2 \ 


n 2 -0 o 2 




0 o +'i 

(45) 


The integral expression for w z=0 in this conical region, both 
on and off the surface, reduces to the same equation which is 


w. 


z =0 


PC Q 0 O p V Mi 2 2 


2rrV n J-x Tix 0 o _1 h 


in: 


(46) 
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Region I 3 : inside both Mach cones and ahead of the plane 


y = - (x-c) . The solution in this region can be produced by 

3 

subtracting from the value of w z _q given for region I 2 the value 
of - w z _q given for region I x except that in the latter case S Q 

is replaced by and x by (x-c). Thus, for this region. 


w 


z=0 


= —° I L j — arc cosh — 

2trY 0 Le o B 


+ •£. 


l-6r? 


arc cosh 


2-9 0 g-T)g 


. 2 

T — 


Si 


cs 

arc cosh 

t ) e x 


arc cosh 


(2-0 o 2 )cu2-v,2 “ 

0 O 2CU2-T|2 


(^7) 


Although the uniformly loaded lifting surface was the only 
prescribed loading analyzed, it should be noted that the basic 
integration leading to a solution of this type of problem (equa- 
tion 0-7)) is in no way restricted to a uniform load. Arbitrary 
loadings that may or may not be analytic functions of x and y 
can be specified and the problem therefore becomes one of 
technique in integration. The solutions for the uniformly loaded 
surfaces, however, are particularly useful. By methods of 
superposition these solutions can be used to obtain the surface 
loading for specified plan forms (the inverse problem) as will be 
illustrated in the following section. 


Load Distributions on FDnt— Plate Lifting Surfaces in 
Supersonic Flow 

Infinite span wing .— Since the vertical induced velocity is 
constant for the supersonic airfoil of infinite aspect ratio (equa- 
tion (22)) and uniform load it follows that the airfoil is a flat 
plate. This property distinguishes the infinite aspect ratio 
problem from all other plan forms considered, for the load 
distribution must be modified in the latter cases so that twist 
and comber are removed from the wing to obtain a flat plate. 

Denoting the angle of attack of the airfoil by a. 


a 


w 


Vo 


2V 0 2 



(43) 
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Moreover, 

Pr-P u - p o (c V^i> " p o c o 

and, setting 

Pl-Pu = 

5-Po^O^ 0 

it follows that 

Ap ^ 2C 0 

*7 ' v 

Eliminating C c between equations (48) and (49) 

Ap 4a 

q Vm--2." 


(49) 


(50) 


The result given in equation (50) is the well-known Ackeret expression 
developed in reference 10. The derivation here follows the approach 
of Prandtl (reference 2). 


Rectangular plan form .— Since the vertical induced velocity 

for the uniformly . loaded supersonic airfoil of rectangular plan 
form is not constant over the portion of the wing covered by the 
tip Mach cones, it is necessary to modify the load distribution 
within this region in order to get a flat plate. The determination 
of the required load distribution will be shorn to depend on the 
solution of an integral equation and subsequent problems dealing 
with other plan forms will, from a mathematical standpoint, be 
similar in form. 

The rectangular plan form will be thought of as being built 
of superimposed trapezoidal lifting surfaces with variable angles 
of rake (fig. 8), each trapezoidal surface having a constant load 
distribution but with loading allowed to vary with the variable 
rake angle 5 . 

Since the flow over the part of the airfoil within the Mach 
cone is conical, it is possible to express w z _q as a function 
of T) where 

= 

x 
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Setting 


0 tan 6=0 


and using equations (42) and (43) 

O' (e) (-«- /'VEE 6.9 ( 51 ) 

0 J 0-0 \ J_ i T)! 'Ql+9-' 

where C'(0) = cp^-cp^ for the single trapezoidal surface with rake 
angle S . 4 

The solution of the problem depends on the determination of a 
function C'(0) which, when substituted in equation (51), will 

yield a constant value of w z _q(t])j that is, a value of w z _q 
independent of the variable tj. Imposing the condition that 

= o 

dTi 

the problem resolves itself into one of solving the equation 

o - i ./' 0 ' 1 cmeT iS? -iii 
to\J'e = o J- 1. % ’ii 

Introducing the notation 


+9 


J. 


r' ] ^ i ~ T ix 2 _ 


Til T il + 9 


&1 (^ , ^ ) 


the integral equation may be written in the form 

0 ■f’Vo r j~[ f>=_ ’ ,_e C(e)5i(n,9)ae + ± P 0=1 c' (elchh.eHel 
Jpi=n -'A=_n + e J 


0= — q+6 

where the singularity in the integrand necessitates the use of the 
infinitesimal € . Taking the derivative, 


0 = 


lim 
e — 5>0 


f 1 C f (0)~ d0+^ 1 c’(0)l^cL0-C , (-n-€)G 1 (Tl,-T}-e) 

L4 » Stj 'J e 9 t) 


5Gi 


+ c’ (-r,+e)Gi(T],-q+6) 
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It is obvious from equations (42) and (43) that, if C'(0) is a 
continuous function. 


lim 
6 — > 


0 


-C'h-£) Gi(n,-4-e)+C’ (-n+e) Gx (T],-ri+e) 


= 0 


Eence 


V 1-11 1 2 /i c 1 (e)'d e 
j 0 nx+0 


(52) 


and the solution can he written in the form 


C'(0) 


•JeJi-e) 


(53) 


where Ci is a constant to he determined later. Substituting from 
equation (53) into equation ( 51 ) 


2rtY o w(r|) z=0 

PCi 


p i ird 0 n i d 6 0 '■*/ i —T l 1 2 

Vs (1—0) Je(i-e) J -i ni(ni+ 0 ) 


dT): 


(54) 


The region of integration in the i'i x , Q plane for the double 
integral of equation (54) is shown as the cross-hatched area of 
figure 9, a singularity in the integrand occurring along the line 
0 = — ‘ Hi. Rewriting the equation and reversing the order of integra- 
tion in the double integral. 


2flV e w(r)) z=0 

PCi 


— 2rt arc 


sin 


Vi 


V_ lim . r r ^h d€) 

Jo C —W _1 T !x '\J 0 (r\i+eyje(i-e) 


r> i d0 

f / _____ 

V_T'4£ (^] X"^0 ) V0 (1 — 0) - 
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By means of the substitution 

9 = **, rji = -p2 

the bracketed expression in this equation can be shown to vanish 
for all values of t ] 1 between zero and —1 so that, finally. 


— CiP r—~ 

w__a = 2 it arc sin v0 

Z 2*V 0 


1 

o' ~2Vo 


(55) 


Since the trapezoidal lifting surfaces are superimposed, the 
loading C(9) over the resultant rectangular plan form satisfies 
the relation 


cL0 sj9 (1—0) 


Imposing the condition that C (6) =0 at 0=0, it follows that 

C(0) =2 Ci arc sin//© - (56) 


This equation gives the incremental change of acceleration potential 
between the upper and lower lifting surface of the rectangular wing. 
As a result . the increment in pressure is 


P = P 
ma o 




2p Ci arc si 
o 


n/lT 


Expressing the pressure difference in nondimens ional terms. 


. £E . i£i arc sto/T (57) 

~P 0 V 0 2 <1 V o 2 

The constant Ci may be eliminated between equations (55) and (57) 
and as a consequence 

Ap W„_A 8 / — 

— £=!£ arc sin-v 9 

<1 v 0 
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Since the angle of attack a of the airfoil is by definition equal 

to the final expression for the loading, in coefficient 

V ° 

form, over the outer portions of the rectangular wing is 


^ = -J*L- arc sin Je ( 58 ) 

a * ,jw~i 

The general approach used to obtain this result is similar to that 
used by Schlichting (ref erence 11). The error in Schlichting’s final 
result has been noted by Busemann and others. 

Lift coefficient C T for an arbitrary wing is defined by the 
relation 


C L 



(59) 


where 


L total lift of the wing 

dS element of area on the wing 


total area of wing 


For the rectangular wing the. values of Ap/q over the tip 
and center sections are given by equations (58) and (50). As a 
result of this integration 


C L = 


Ua 

P 


1 

2 pA 


(6o) 


where A is the aspect ratio and by definition equal to the ratio of 
the square of the span and the wing area. As a final conclusion the 
lift curve slope of the wing is 


dO 


L _ 


da 


P 


1 - ■— ) 
2pA/ 


(61) 


Trapezoidal plan form .- The results given in equations ( 58 ) and 
( 6 l) are capable of generalization to the case of the flat plate 
having trapezoidal plan form and with rake angle Oq less than the 
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Mach angle of the stream. For such a configuration the airfoil is 
again blanketed in part by the tip Mach cones and the loading in 
this outer section of the airfoil must be adjusted properly to give 
constant induced vertical velocity. Superposition of trapezoidal 
lifting surfaces with loadings varying with rake angle 5 can 
again be used and the conical nature of the flow employed. Setting 

n = Py/x 

0 = p tan 5 

6 n = p. tan F> 
o o 

equations (42) and (1+3) give the expression 


2irV o w(n) z=0 

- 



v J-1 'll 


r iX +eS 


de 


(62) 


where C f (0) = cp^-®^ the single trapezoidal surface with rake 
angle S. 

The analysis in this case follows along lines directly analogous 
to that used for the rectangular surface. For the present configura- 
tion the loading function for the superimposed trapezoids is given by 
the relation 


C»(0) = ^ ■ - 

J (Q—6 0 ) ( 1 — 0 ) 

and by means of the substitution 


(63) 




1—0 o 


the integration to obtain w can be simplified to give, as a 
final result, Z "° 


PCi p l ■ nd0 _ 3Cx« 

2nV 0 J. /”(0-0 o ) (1-07 2 Vq 


m 



MCA TN No. 1412 


33 


The loading C (9) over the resultant trapezoidal plan form 
can he found from the relation 

dC(0)_ Ci 

~Ie~ ~ J {e-e 0 ) (i-4) 

From the "boundary condition that C(d) = 0 at Q - 0 Q it follows 
that 


and 


c(e) 


2Ci arc 



Pl-Pu 

fPoV 


q 


4Ci 

v o 2 


arc 



(65) 


Elimination of C x between equations (64) and ( 65 ) and introduc- 
tion of angle of attack a for -w z=0 /V o gives as aerodynamic 

loading over the portion of the airfoil within the tip Mach cones 
the expression 


Ap 8a 

— = — arc sm 
q 



0<e o < 1 


(66) 


Figure 10 indicates the variation of the loading over the tip section 
of the trapezoid. The variable (p/a) (Ap /a) is plotted against 
p tan 5 for P tan 8 0 equal to 0, 0.3, and 0 . 6 . The curve for 

P tan 5 0 = 0 corresponds to the case of the rectangular' wing and 
shows results in agreement with equation ( 58 ) . 

By means of equation (59) together with equations ( 66 ) and (50) 
the lift" coefficient of the trapezoidal wing is expressible in the 
form 


c L = 


4a 

T 


/i- 


c 

2b 


tan 5 r 


tan 

2b 


h\ 


7 


1 — £• tan 6 0 


( 67 ) 
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Introducing the aspect ratio A of the wing where 


c (1 - £ tan 5 0 ) 

lift coefficient can also be written 



From equation (67), 

dC L U A ~ tan 5 o ~ |b tan M 

da 3 \ 1 - ~ tan b n 

D u 

In figure 11, 0 . is plotted as a function of A0 for 9 = 0, i, 

act <= 

and 1. The curve for 9 0 = 0 agrees with results given by equation 

(6l) for the rectangular wing. All curves are terminated at values 

of A0 for which the tip Mach cones intersect on the trailing edge 

of the wing. 

Triangular plan form, type 1.— Hie pressure distribution over 
triangular lifting surfaces with constant induced vertical velocities 
will be developed in the following three sections. These plan forms 
are indicated in figures 12(a), 12(b), and 12(c) and shall be denoted 
respectively as types 1, 2, and 3. Types 1 and 2 are actually special 
cases of type 3; namely, the cases where one leading edge is parallel 
to the free stream, and where both leading edges make equal angles 
with the stream direction. Type 3 includes any plan form which has 
leading edges swept behind the Mach cone but on opposite sides of 
an axis drawn through the vertex of the triangle and parallel to 
the free stream; and, further, has a trailing edge such that the 
Mach cones from either tip do not cross the surface of the wing. 

The principal reason for considering the three types separately 
is to show the manner in which the spanwise loading appears in 
the solution of the problem. In types 1 and 2 the proper load 
distribution is found readily while the final type requires a more 
careful treatment. 

In order to determine the load distribution over the airfoil 
it will be convenient to use a differential element over which the 
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loading is constant. T!ie elements may then he summed and the distri- 
bution of loading adjusted so that the induced vertical velocity at 
any point on the total lifting surface is constant. For the tri- 
angular plan forms it is possible to assume that conical flow exists 
and the analysis may be carried out using the angular coordinates 
which have already been introduced. 

Figure 13 shows the elemental lifting surface to be used. The 
sides of the element extend back from the tip of the Mach cone; 
making angles & and & + AS with the positive X-axi3 or free— 
stream direction. Corresponding to previous notation, the relations 
9 = 3 tan S and 9 + A8 = 0 tan (& + AS) are used. The vertical 
velocity induced by the element of surface may be denoted by Aw and 
it follows that 

Aw = w (3+A0,T)) — w(9,tj) 


where w(8,i]) and w(0+AS,tj) are the velocities induced by the 
triangular— tip surfaces with constant loading and with tip angles 
equal to S and S + AS, respectively. Applying a limiting 
process. 


lim 


A9 — ^0 


Aw 

A9 


w(6+A0,ti)-v(0,t)) 
lim 



&w 

be 


(70) 


It follows that w z _o for the resultant lifting surface will be 
evaluated by an integration with respect to 9. If bw z _q/^ 9 can 

be expressed in the form of an integral with respect to r^, the 
relation for w z _ 0 will then be similar to those given in equations 
( 51 ) and ( 62 ) for the previous plan forms and the expectation will 
be that the function C (9) can be determined to give constant 
induced vertical velocity. 

The method of attack just outlined is postulated on the 
existence of an integral expression for <3 w 2 _q/< 30. Such an 
expression is obtainable directly. from the integrals in equa- 
tions (33 ) , ( 36 ), and (39); however, the alternate relations 
given in these equations lead more directly to the proper form. 
Proceeding to this calculation, the following results are obtained 
for the elemental lifting surface; 
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For -0< i) <0 


^2=0 = 

30 2nV 0 


1 -11 
— arc cosh 

e 2 


arc cosh 


^ 02 a/1-02 


1 + 0 T | - v / l - TJ 2 " 
0+n . 0(0+T)) . 


( 71 ) 


for — 1 < Tj < — 0 


dv z=0 PC! 


30 2itV r 


1 -1 

— arc cosh — 

L0 2 tj 


02 71-02 


1 -(1+011) ' 
arc cosh - 


0+Tj 0 (0+TJ ) 


(72) 


and for 0< tj < 1 


3v 0 pCi I i , i i v 1+0T i Vi-n 2 

= — arc cosh arc cosh 

30 2jtV 0 L 0 2 


n e 2 a/i— 0 2 


0 +n 0(0+n) J 


( 73 ) 


It can he verified easily that 


3 . -1 -1 

— arc cosh — = — — • = 

3n n tj 7 1 ~ T 1 2 


for tj < 0 


and 


while 


3 1-1 

— arc cosh - = — f== for ^ >° 
3rj ■= il ii Vl-i] 2 


3 " . 1 + 0 T ] 1 / 1—02 

- arc coah - —J — 


for -e <>i <1 
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and 


Sr, 


arc 


(l+9r,) 1 /II02 . . „ 

cosh -1 = / for — 1< p , 

e+M e+r> V 1-TJ2 


Combining these results, the desired expressions result: 
For T) < 0 


Sw z=0 PCj. f - a/1 -n 2 ^ 3 f -1 


* 


2«V 0 1 6 (0+H ) sLi '-0 2 '1ia/i-'H 2 <9 2: (s+n a ) 


dr.] 


-PCi 

' Vl-n 2 

d'Qn 

2rtV 0 

Le( 0 +n) J_ x i 

?»*{ 1 i — Tj -yp _! 


(74) 


and for tj > 0 


dw z =Q = — PCi f Vl-jg + r n 
oa 2 rf 7 „Le(e+r‘ 


d-ni 


0 La(a+Ti) eii i (a+'ii) 


'll 2 j 


(75) 


Summing the elements over the type 1 triangular wing, induced 
vertical velocity is 


w ( t) ) 0 = _±_ r °c(e) r + P S!Ll " 

2n ^J 0 L 0 ^ +r l) 0 tj 1 (0+^ 1 ) a / 1-^ 1 2 - 


d£> (76) 


where C(0) = cp u -<p^ for the element at G = arc tan — , If v s _q 
is constant, then 




z=0 


= 0 


dii 


and from this criterion the function C (d) will be determined. 
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Thus 


_6 r 8 r V l-u 2 C(0) d 
o = _JL — / : — — de + 

2rtV 0 L dr\Jo 0(9+ tj) 5' 


2L_1 

iW 0 0r)i (e+.ii) ,j l- 1 ] l 2 - 1 


and, after introduction of the function G2 (h,0) where 


G2(tj,0) » — ■ . ^ 1. — 

l 0il 1 (0+ r l 1 ) 1 


the integral equation may he bitten in the form 


0 = vl— tp / 


; /->$>C(0)d9 


St] 0(0+T|) 


+ lim 

e — >0 


[rcw^^fcw^d, 1 

Wo ^ T > '-Ln+e ^ T 1 


00 


1 -1 m P *" 

+ ^ ^ ^ C ( — 1— € ) G2 (’: , — 'l — £ ) + C(— T)+e) Gp (1] , 11+6 ) 


(77) 


The function Ga(i),0) is given hy the first two terms in the right- 
hand members of equations (71) and (72) and from consideration of 
these terms it follows that, provided C(0) is a continuous function, 
equation (77) becomes 


0 _a./T^r 9 o c«m. 1 r b m de (7 e 

5 ti J 0 (0+rj) fj N 1-;] 2 J 0 0(0+11) 

Comparison of equations (52) and (78) , together with equation (53 ) , 
indicates that the function C(0) is of the form 


C ( 6 ) 


Cx 



(79) 
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where Cx is a constant. Substituting from equation (79) into 
equation (76) . 


e o cL0 T 11 dxii 

Ve(e o -e)^-i ix aWW 

(80) 

The region of integration in the i]x,0 plane for the double 
integral of equation (80) is shown in figure 14 for the case in 
which — 0 o <ti<O. A singularity in the integrand of the double 
integral exists along the line 0 + 'Ox = 0. Reversing the order of 
integration in the double integral, equation (80) may be rewritten 
as 


2jtV w 
o 

pc 


Vo = P 

1 J n 


Vl-T| 2 d0 


o (0+h) J 0(0 O - 


= -/ 

I— 0T J o 


2itV-W__ n 1 p ®o 

PCx J ° 


d0 


(e+n) J o(V 0) 


i 


" e o drix p®0 


d0 


X i'll «/ l-Tix^o (0+Tix) a/ e(0 o -0) rix Vl-TJx 2 ^o (0 +Tx) V 0(e o -0) 


r 

^ r» _ */l -n _ 2 'JQ 


6o 


d0 


(81) 

The single integral in equation (8l) has a singularity at 0 = — t) 
since — 0 o <r,<O and tj therefore lies inside the region of 
integration. A corresponding singularity occurs in the second of 
the double integrals at 0 = — Tx. Consider, therefore, the integral 


l 


9o lim 

o (0+1]) V0(0 o -0~) e — > 0 




d0 


d0 


u w o 


___.+ r w o _ 

(0+T1 ) P~Q 0 -e) J ^ (0+T1 ) P 0 (0 o -0) J 


-“T+G 


The indefinite integral is 


(82) 


J — Tj0 — 1 r\ 2 
o 


In 


-n0 o +00 o +2n0 - 2 V(~T]0 o -T} g ) 0 (0^0) ' 


6 +tj 
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so that, the definite integral is 


lim 


1 ©2 ^ Tll [- 2ne Q-€ a o- 2T l 2 - 2eT l- £ ' v/ (-^Q-n g ) M-e ) (3 0 +T i + g) 1 

0 jzppp? | “®o [-2T>e 0 ^ee-2ri 2 +2€Tr : 2 A / (-ns 0 -n 2 ) (-11 +eX0 o +r i-«O 1 


The value of this expression is 0 and equation (8l) can therefore be 
written as 


2ttY o v z---0 

3C X 


h 9 ° dT H f 

■J_ x n,i v l-pjo 


e 

0 


d0 

(e+rix) Ve(0 o -e) 


(83) 


Since, in this region of integration, — l<i)<— 0 O it follows that 



d9 

(0+il 1) J 0(0^0) 


1 

■ — , arc 
^1^0 +T, >1 2 


tan ~ ril9 o +0e o +2 r ‘ l9 
2»/ ( t !i 2+1 1i0o) ® (®o~ 



—IT 

*/i(e 0 +fi) 


and 


w 


z=0 


3Cx P~® °' 1 

OT o li ’ll -A-iiV 'ii(ii + 8 0 ) 


(8k) 


The integral of equation (84) can be transformed by means of 
of purely algebraic substitutions into a form which integrates 
immediately into complete elliptic integrals of the first and second 
kind. However, in the consideration of the type 3 plan form it will 
be necessary to resort to other methods of transformation, so that a 
more uniform approach, employing Jacobian elliptic functions, will 
be used throughout. Reference 12 contains a complete discussion of 
these functions. 
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The quart ic under the radical in equation (84) can he. reduced 
to an expression of the type appearing in elliptic integrals of 
canonical form. This is accomplished by successive application of 
the transformations 


= - 


k+lS 

1+S 


and S 


k 

t 


where k and l are chosen so as to destroy the odd powers of the 
variable. By means of these transformations, induced velocity 
becomes 


(1— k 2 )^ 2 ok + (k— ) t— 1] d o 

27 o VF(e 0 -k) (1— k 2 ) (1— k 2 t 2 ) a/ ( 1— t 2 ) (k 2 t 2 — l) 


(85) 


where 


1 _a/i_ 0 o 2 


The integration of equation ( 85 ) will be performed after first 
considering two parts such that w 2 _q = Vx +w 2 where 

_ pCx (l-k^k 2 /'' I (k--)tdt 

2V o Vk(S 0 -k) (1-k 2 ) J 1 (i-k 2 t 2 ) J (l-t 2 ) (kH 2 — l) 

and 

1. 

■PC-1 aEI (1— t 2 )dt 

W2 — — ■ ■ . ■ ■ ■ — / — 1 , ■ ■, — . r _- 

2V o V k ( 0 0 — k ) ( 1-k 2 ) J x ( 1-k 2 1 2 ) a/ ( 1— t 2 ) (k 2 t 2 -l ) 

This separation is prompted by the fact that the integral for 
w x is expressible in terms of elementary functions after the simple 
transformation t 2 =z. The results of such an integration lead to a 
value which is zero at the lower limit and infinite at the upper limit. 



42 


NACA TN No. 1412 


However, an inspection of the original integral in equation (84) shows 
that w 2 _q is finite so the infinity obtained for w x must he 
canceled by a corresponding infinity of equal magnitude in w 2 . The 
actual proof of this statement necessitates, of course, treating the 
combined expressions as an indeterminate form where the upper limits 
of the integrals for w x and . w 2 are replaced by A + e and the 
limit is taken as e approaches zero. £ 

Introduce now in the integration of v 2 Jacobian elliptic 
functions and set 


t = sn (u,k) = sn u 


so that 


dt = cn u dn u du 
The expression for w 2 becomes 


v 2 


PC X (l-k 2 )k 2 i ynK+iK* 

2V o ^k(e 0 -k)(l-k 2 )Jk 


cn 2 u 

dn 2 u 


du 


where K and K 1 are the complete ellip tic integrals of the first 
kind with respective moduli k and k!«vl— • k 2 . Integrating and 
combining with w x , it follows that 


w z=0" iW i + 


P£i (l-k^k 2 ! JL_ 
2V 0 V k(9 0 -k) (1— k 2 )Lk 2 


u + k 2 


sn u cn u 
dn u 


s-j K+iK’ 

/Ik 


where E(u) is the incomplete elliptic integral of the second kind. 
After . substitution of the limits. 
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where E* is the complete elliptic integral of the second kind with 

modulus k’ = Vl-k 2 . Equation (86) can be further simplified by 
writing k in terms of 0 Q so that 


w, 


z=0 


/ 

= _ pci j 2I1+ 1KI 


2V, 


5 o 2 


E' 


where the modulus of E‘ is ^/l-k 2 and k = 1 ~ .^^o 2 

G o 

For the loading in question 

p r p u ■ 'VvV ’ <>o 0 vA" 


v ©'q-o 


( 87 ) 


and 

Pl"Pu = 4P _ 2Ci / Q 

|Po V o 2 " a - ” W ^ 


By mean .3 of equation ( 87 ) the constant Cx 
since a = - w 2= 0 /V o , 


Ap 

<1 


2a ke (1- 
be 1 J e 0 -e 


may be eliminated and. 


( 88 ) 


Figure 15 shows the variation of E with 

a q 

3 tan & Q equal to 0.?, 0.6, and 0.9. 


3 tan & for values of 


From equations (88) and (59) the lift . coefficient, of a right 
triangle wing can be determined. It follows that 
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Since the aspect ratio A is given hy the equation A = 20 o /3, 

equation (89) can he used to find the lift curve slope as a function 

dCr 

of A. In figure l6 a plot of 3 as a function of PA is given. 

Triangular plan form, type 2.- Figure 12(h) shows the symmetri- 
cal type of triangular plan form considered in this section. The semi-r 
vertex angle is S 0 and 0 Q is defined hy the relation 


0 O = p tan & Q 

The loading element used in the previous section can he used again - 
and equations (74) and (75) are applicable directly. Because of the 
symmetry of the figure, it is necessary merely to insure the constancy 
°f w _ 0 over the left half of the wing in order that the entire wing 
he specified a flat plate. 

When the lifting elements are on the left portion of the figure, 
the induced vertical velocity contributed hy them will he given hy 
equation (76) for, aside from the change in the loading function 
C(0), this situation corresponds exactly to the condition considered 
in the type 1 plan form. To equation (7 6) must also he added the 
expression for the induced vertical velocity on the left side of 
the wing produced hy the lifting elements on the right. In the 
development of equations (74) and (75) the element was situated in 
the second quadrant of the axial system; that is, y was negative 
ap ri x was positive. For an element on the right side of the wing, 
equations (74) and (75) are again applicable if the axial system is 
changed so that the Y-axis reverses its direction while the X-axis 
remains fixed. Since, hy definition, 

TJ = P^ 

1 y X 


it follows that the change of sign for y results in a change of 
sign for TJ. When the lifting element is on the right-hand side of 
the wing the egressions for gradient, of induced vertical velocity 
are, thus. 


for b <0 


,^fz=0 _ 

be 


PCi 

r a/ i-ti 2 x 

pT\ diu 

2*V 

- e(e-ri) J. 

_ x ' 'll) aA -t 1i 2 - 


( 90 ) 
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and for t]> 0 


^ v z=0 

ae 


pCi r a/i— , n n dT| i 

2«v 0 L 8 (0— tj ) J i 8nx(e-n x ) - 


(91) 


Summing the elements over the type 2 triangular wing, induced 
vertical velocity over the portion of the wing for which -1< < 0 

is 



f vi-n 2 

[_0(6+n) 


► r 

j 


n 


dtji 

©hi (8+%) V 1— r i l 2 


C(0)d0 


g p 9 * r + r n am 

2jrV oJo L 6 J_x 0T, x (9-T) x ) Vl-h! 2 


C(9)d0 


2nV 0 ^ 


r 


dT]] 


a/1-v. 


L02 -“' r > 2 -i m(e 2 -m 2 ) Vi-nJ 


2 J 


C(0)d0 (92) 


The function C(0) in equation (92) must give a cons tent value for 

w(n) „ so that Sw „/an will vanish and also he such a function 
' 1 ' z=0 z=0' 

that the loading will he symmetrical about the line 0=0. Imposing 
these conditions it can he shown that a solution is given hy the 
relation 


C(0) = (93) 

V0 q 2-0S 

and, after substitution in equation (92), 


ffY o w z=o = r e ° £ 1-P 2 d9 r°° d9 dr ii 

PCx J 0 (0 2 --ii 2 ) V 0 o 2 -0 2 i V 0 O 2 -0 B ^L i Ox Vl-fi 1 2 (02_q 1 2) 


(94) 
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The integration of equation (94) is to be performed under the 
assumption that — 0 O < t)<0 so that the region of integration in 
the 1^,0 plane for the double integral is as shown in figure 14. 

A singularity of the integrand exists along the line ' G + 'Ai=0. 
Reversing the order of integration in the double integral, equation 
(94) may be written in the form 


, _ vC?/' 9 ° as 

P C 1 (02_t,2)a/0 o 2-02 


(95) 


n- 6 ^ dr] 1 n 

j-l % Vl-r a 2 ^0 


d0 


_ r" 1 dT ii r B ° 

J- d o n i J i-0i 2 ( e 2 - Ti 


d0 


fj 1 — T| o (02_rj^2 ^ a/ 0 q 2 —02 tj^a/ 1— rj x 2 o ( 9 2 —t\ yj 0 q 2 — 0 2 


Evaluation of integrals of the form 




0o 


d0 

(02_rj2) J q^Z_qZ 


can be accomplished by means of the substitution 


x = 


e 


After substitution, the integral becomes 


I = 


hf 

3 o 2 Jo 


dx 


On 2 Jo f Jlf, 

? o 2 

and, by straightforward integration, 


(~- - x 2 ) Vl-x 2 
9 0 2 


0 for -! — < 1 

0 2 


e o 2 


>1 


2f\ J n 2 -e 0 2 


for 
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This result shows that the second double integral of equation (95) 
vanishes (since © Q 2 >r i;L 2 ) as does also the single integral in the 
equation. For the remaining double integral, however, Q q 2 < t )i 2 
and 


v zp0 . fii r B o 


dt]] 


j 

~ Oo _i o i 2 a/ i— n i 2 aAi 2 -9 0 £ 


(96) 


Setting 


equation (96) transforms to 


z = — 


JL 


w 


z=0 


= _ £Ci r(?o 


z 2 4z 


2V 0 4 J ( Z 2-l ) (l-©0 2 z 2 ) 


Introducing the modulus k = 0 O and making the substitution 


z = sn (u,k) = sn u 
the expression for v z= q becomes 


w 


z=0 


ipci p 

2V./, 


K+iK' 


sn 2 u du 


K 


iPCx . . . K+iK' 

r [u — E(u)]„ 

2V„k 2 * 


ftCi 

2V 0 k 2 


E 1 


where the prime again refers to the complentary modulus k* -J 1— k 
of the complete elliptic integral. Since k = 0 Q 

PC! 


w 


z=0 


^o 6 i 


• E 1 


(97) 
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where k f = aA-Oq^ 


For the loading in question 


P2-Pu _ AP _ 

l PoVj °^ <i " v 0 2 Ve 0 2 -e 2 


so that, eliminating Ci between this equation and equation (97 ) > 

Ap 


Uo0 o ? 


(98) 


<i jwe 0 2 _e 2 e‘ 

Figure 17 shows the variation of - — with 0 tan & for values of 

a q 

0 tan 5 0 equal to 0.3, 0.6, and 0.9. 

Using equations (98) and (59) it is possible to find the 
expression for lift coefficient of a triangular or delta wing. Thus: 


°L = 


2rt0 o a 

0E* 


Since aspect ratio of the wing is 


(99) 


A.J& 

P 


lift coefficient becomes 


r _ *ccA. 
C L " ?F 


where the modulus of E’ is k* 




A 2 0 2 

IT- 


(100) 

This result agrees with 


that obtained in another manner by Stewart (reference 8). In figure l8 

dC L 


a plot is given of 0 


da 


as a function of 0A. 


Triangular plan form, type 3 . Figure 12(c) shows the plan 

form now to be considered. Belative to the X— axis or free— stream 
direction the sides of the triangle form the angles 6 0 and 61 so 
that the total vertex angle is 5 0 +5i=2A. The variables 6 0 and 6x 
are also introduced satisfying the relations 0 O =0 tan 5 0 > 

e^p tan 61. The same loading element which was used for type 1 
and type 2 triangles may be used and equations (74) and (75) a PPly • 

It will then be necessary to determine the distribution of load so 
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that the induced vertical velocity over the plan form is a constant. 
Since this induced velocity must he the same on both sides of the 
5=0 axis, two equations result: 

for — 1 < r|< 0 


p r 9 ° 

2jt7’ 0 


' V I-’} 2 

l -0(e+q) 



dr) x 

^i(0+'U) 



c(e)de 


( 101 ) 


for 0<rj<l 



- /\/l— p ( 'l 

-6(e+q) J x 


&)! " 

0ni(0+ni) Vi-\ 2 - 


C(9)d0 


(102) 


From the solutions to the problems of type 1 and type 2 it is 
possible to construct a solution of the more general problem by 
expressing the loading function in the form 


c(e) = 


AO + B 

\/~(0l+G) (9 0 -9) 


(103) 


where A and B are constants which can be determined in terms of 
w z=0 f^om equations (101) and (102). Substituting from equation 

(103) into equation (101) and (102), the expressions for induced 
velocity become 


w z=0 = ~~ [AH X ( 00 , 03 .) + BH 2 (0 O ,0 1 )3 


vhere 


w z=0 ~ ^r-( _AH i( 0 ii 9 o) + BH 2 ( 9 i, 0 o )] 
^ v o 


Hi(0x,0 o ) =P 61 ■ Ani 

A ’ll V (l— ' 'll 2 ) (fix — 9 i) ( ,f ii +9 o) 


(104) 


(105) 
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dJJi 

’ll 2 *J 


(106) 


The evaluation of Hi (0i,0o) and H 2 (0i,0o) is accomplished in the 

same manner as has "been used previously: first, a reduction of the 

quartic under the radical to canonical form and second, transforma- 
tion "by means of Jacobian elliptic functions followed by direct 
integration. Since the calculations for both equations are quite 
similar, only in the case of Hi(0i,0 o ) will the details be mentioned. 
All relevant information can be found in reference 12. 

Applying the transformations rji = ~ ~ and S = - , where 

1+S t 


l—0o^i (1— 0o 2 ) (l— ®i 2 ) 

0 O + 0i 

1-0001 W (1-0Q 2 ) (1-0^) 

e 0 + 


(107) 


(108) 


and introducing the symbols 5 and k defined as 

(b— a ) 

H = — • 

V (l-^ 2 ) (0 o -a) (0i+a) 

' k = h£^ 

0 o -b 

equation (105) reduces to 


(109) 


(no) 


1c [t 2 +t(a-b)-]j dt 

(1— a 2 t2 ) J ( t 2 — 1 ) (1— k 2 t 2 ) 

The integrand divides naturally into two parts, one containing even 
powers and one containing odd powers of t in the numerator. The 
latter part integrates into elementary functions after substituting 
t = u 2 and equation (ill) thereby becomes 


Hi = 


a 2 Bk 
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Hi = a 2 Rk 


L k^-a 2 ^ ( 


(t 2 — l)dt 


(l-a 2 t 2 ) JW-L) ( 1— k 2 t 2 ) 


( 112 ) 


Setting 


I2 



(t 2 — l)dt 

(l-a 2 t 2 ) */ (t 2 — 1) (1— k^t 2 ) 


and substituting 


it follows that 


x = an (u,k) = sn u 


I2 


-i/ 

J K 


K+'iX*' 


1 + (a 2 -l) 


3 n 2 u 

l-a £ sn 2 u- 


du 


If sn r s r , Ig now may he written as 

dn 7 


I® « -K* - i 


k 2 sn 7 cn 7 K 


pL iA £? 3 n y: C n 7 in 7 sn 2 u 

/ — du 

./t? 


1— k 2 an 2 7 sn 2 u 


or 


Is = -K’ - 


V 


A 2 - 1 r„ 




ii (u, 7 ) 


-j K+iK* 
_ K 


(113) 


where II (u, 7) is the fundamental elliptic integral of the third kind. 

The evaluation of H(u, 7 ) is best achieved by means of its 
expression in terms of theta- functions and zeta functions. Thus 


- n (u,y) - £ log + uz(y) . 

2 © (u+r) 

and, substituting the limits in equation ( 113 ) , 


n(K+iK‘,7)-4I(K,7) 


1 log ©(K+IK*— 7) Q (K+ 7 ) + 

2 ©(K-a) ©(K+iE’+r) 


iK'Z (a) 
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The Theta -functions are a class known as quasi— periodic, that is, 
they satisfy the relations 


0 (u+ 2 K) = 0 (u) 


-(i:’-iu) 

0(u+2iK’) = — e K . 0(u) 

From this property, together with the fact that ©(u) is an even 
function, it follows that 


1 ft? 

6 (K+iK*-y) © (K+y) _ k 
©(K-y) 0 (K+iK’+y) 


Moreover, since 

z( 7 ) = E(y) - r | 
is. 

there results 


n(K+iK',r) -H(K,7) = IK’Efr) - i7 ~+i ~ 


and 


I 2 = — K* + { — + K« [e(7) - 7 -1 } (114) 

v k 2 -#. 2 L2K L KJ J 


The expression for equation (112) can now "be written 


Hi = a 2 Rk [till f 1 
* k 2 -a 2 \ 2 


I - 1 

K 


+ K» 


E(y)— y - 
KJ 


} 


-a 2 RkK* 


(115) 


where the moduli are k for the nonprimed. functions and k * 


for the primed functions. By definition, 7 = arc sn ^ = K (pk) 
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where F i3 the incomplete elliptic integral of the first kind with 
argument — and modulus k. 

1C 

Using the same notation, the equation for H 2 is as follows : 


H 2 = 


KRb 

b 2 }^. 


~{a 2 (l_k 2 )K’-(l-^2)E»-a(l-k 2 ) /^|K'E(r)- 7 -' +- (- - l)l [ 

-ll V h 2 k 2 — 1 l K 2 NEC - J J J 


( 116 ) 


Formulas (115) and (116) can now be combined with equations (104) 
to give 


and 


where 


A • - T5U 2 J\ 


_2G_ 

* 9 O +0 1 


b=- Ip - 1 z - . P , 2e o 0i / - ■ —■ ■■ 

0E* V 0 O +01 


G » 


1 + 0 O 0! - /7l-0 o 2 ) (1— ©I 2 ) 


0Q + ^1 


(H7) 


(118) 


(119) 


and E* is the complete elliptic integral of the second kind with 
modulus V 1-G S . 

From equations (103), (H7), and (118) 


Pl-Pu 

§P 0 Vo 2 


£E = gC(0) = _2_ A9+B = 2a_ / 2G (90-01)9+2000!) 
q T 0 2 Y 0 2 5/(0!+ 9) (0 O — 0) PE* J 0 o +0i V (0i+0) (0 O — 0) 


( 120 ) 



It should "be remarked that the slope of the loading curve is zero at 

e *» 0. Figures 19(a), 19(h), 19(c) show the variation of - — 

a q. 

with 3 tan 5 for values of 3 tan &i equal to 0.3, 0.6, and 0.9, 
respectively and for 3 tan & 0 equal to 0, 0.3* 0.6, and 0.9. 

From equations (59) and (119) the lift coefficient for a type 3 
plan form is obtainable. Two cases- will be developed here: first, 
when the trailing edge of the wing is perpendicular to the stream 
direction; second, when the trailing edge of the wing is perpendicular 
to the line of symmetry. The first configuration may be referred to 
as a skewed wing while the second configuration may be referred to 
as a symmetrical delta wing at an angle of sideslip. There results 
for the skewed wing then 

Cr = — / (e o +0 i)2G (121) 

u E'3^ 

where G is given by equation (119) and E' has the modulus V 1— G 2 . 
This result agrees with that given by E. C. Eoberts in an abstract 
in reference 13 . 

For the more practical case of the delta wing at an angle of 
sideslip, figure 12(c), the lift coefficient can be expressed as 


n 2cut * 

C T - sec A 

L E* 


tan A 


■3 


( 122 ) 


where A is the angle of sideslip and PA the angle between the 
leading edges, E' still has the modulus Vl-G 2 , and G is' expressed 
in terms of 9 0 and. 9 \ which are, in turn expressed in terms of A 
and A by the following equations 


0 O « 3 tan (A4A)' 
6>i = 3 tan (A-A)- 


(123) 


It should be noted that since the pressure distribution has been 
computed only for wings -with.: leading edges behind the Mach cone 
springing from the apex and with a trailing edge ahead of the Mach 
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cones springing from the wing tips, formula (122) is valid only for- 
cases where 

H + A < 90 ? 

A + A < |a [ (124) 

A - A > 0 J 


These restrictions are practically always met, however, for angles 
of sideslip which are likely to be encountered in flight. 

Equation (122) is plotted in figure 20 where 0 i s shown as 

da 

a function of sideslip and A* The figure shows that up to I 5 0 of 
dCr, 

sideslip 0 remains practically constant, 

Ames Aeronautical Laboratory ^ 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 
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Fig. 2a 
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(a) Rectangular plan form. 

Figure 2.- Mach forecone from point P(X,Y,Z) intersecting surface 
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Fig. 2b 



(b) Triangular plan form. 
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Fig. 4 



Figure 4.- Lifting surface with Mach cones and coordinate system for rectangular plan form. 
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